import scipy.stats as stats
import numpy as np

p = 0.3 #客户购买概率
n = 6 #客户数
X = np.arange(0,n+1,1)
plist = stats.binom.pmf(X, n, p)
print('恰有4个客户购买的概率:',plist[4]) #恰有4个客户购买
print('购买的客户不超过一半的概率:',1-(plist[4]+plist[5]+plist[6])) #购买的客户不超过一半
print('至少有一个客户购买的概率:',1-plist[0]) #至少有一个客户购买